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Summary 

Analyses  are  made  of  theoretical  seismograms  of  torsional  oscillation 
caused  by  a  localized  surface  source  for  a  homogeneous  sphere  and  a 
homogeneous  mantle  with  a  liquid  core.  Special  care  is  paid  to  the 
effect  of  parameters  cm  the  accuracy  of  analysis  and  the  limitation  of 
the  applicability  of  various  methods.  At  first,  a  simple  description  is 
given  about  the  assumptions  under  which  the  theoretical  seismograms 
were  compute<^8setion  i).  Next,  the  spectrum  analysis  of  the  seismo¬ 
grams  is  carried  on  by  the  usual  method  of  Fourier  transform,  giving 
useful  information  about  how  to  employ  the  method.  This  analysis  corro¬ 
borates  recent  similar  analyses  of  the  Chilean  shock  of  May  196^, 
(Session  a).  Group  and  phase  velocities  calculated  by  methods  ordinarily 
employed  show  good  agreement  with  theory.  Calculations  of  initial  phase 
also  agree  well  with  theory,  giving  confidence  to  similar  studies  using 
waves  from  earthquakes  and  cxplosion^fisotion  3).  The  idea  of  an 
equivalent  surface  source  is  proposed  relative  to  the  problem  of  the  depth 
of  the  source  within  the  Eart^Seetmt-4). 

1.  Introduction 

In  a  previous  study  (Satd,  Usami  &  Ewing  1962)  disturbances  were  numerically 
calculated  at  various  points  cm  the  surface  of  a  sphere,  assuming  a  torsional  force 
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on  the  surface.  Since  this  it  a  kind  of  theoretical  seismogram,  it  can  be  assumed 
to  be  relatively  free  from  noise  and  from  uncertainty  related  to  the  source  mecha¬ 
nism  and  the  dispersion  relation.  Consequently,  by  treating  this  seismogram  as 
if  it  were  an  ordinary  one  observationally  obtained,  it  will  be  possible  to  find  where 
errors  are  introduced  by  the  standard  techniques  of  analysts  in  common  use,  to 
determine  the  precision  of  each  step,  to  measure  the  effect  of  various  parameters 
on  the  accuracy  of  analysis,  and  to  estimate  the  limits  of  the  applicability  of  these 
analytical  methods.  The  results  of  these  analyses  are  closely  related  to  several 
recent  works  relating  to  the  Chilean  earthquake,  as  well  as  to  studies  of  the  data 
from  earthquakes  and  explosions. 

The  theoretical  seismograms  used  in  this  paper  were  calculated  under  the 
following  aasumptions: 

(1)  The  medium  is  either  a  homogeneous  elastic  sphere,  or  a  homogeneous 
elastic  mantle  with  a  liquid  core  (Satd,  Uaami,  Landiaman  &  Ewing  1963). 

(2)  Uniform  shear  stress  is  applied  on  the  surface  in  such  a  way  that 


8<0 1,  6>$t 

O\<0<02. 


(3)  The  time  dependence  of  the  source  function  is 


(*•*) 


t<h,  t>h 


(1.2) 


(4)  Axial  symmetry  around  the  pole  is  assumed,  hence  m  =  o  and  only  the 
azimuthal  component  exists. 

Calculations  are  carried  out  using  non-dimensional  quantities,  but  for  clarity 
and  convenience,  the  most  frequently  used  quantities  are  tabulated  in  terms  of 
ordinary  units: 

a  (radius  of  the  sphere)  =  6370  km 
2 ira  (circumference  of  the  sphere)  =  40000  km 
V,  (velocity  of  S  wave)  =•  6-667  km/s. 

As  a  result  of  employing  the  above  values,  it  is  easily  seen  that 
(2 ira)IV$  (unit  of  time)  =  100  min 
Duration  of  the  shear  stress  being  applied  =  3  min 
Diameter  of  the  area  where  the  stress  applied  —  444  km. 

The  following  notations  are  used  in  this  paper: 

iTu :  Torsional  mode  of  the  free  oscillation,  i  and  a  being  the  radial  mode 

number  and  the  order  of  the  oscillation,  respectively. 

(ton:  Azimuthal  displacement  due  to  tTn  mode 

(<»«):  $  ?  <w« 

/ (t):  Time  function  of  the  force  applied 

f*  (p):  Fourier  conjugate  of  the  function / (t). 


2.  Spectrum  analysis 

The  theoretical  seismogram  is  the  superposition  of  the  various  free  oscillations 
having  incommenaurable  discrete  frequencies,  while  the  analysis  is  based  on  the 
Fourier  transform  of  a  continuous  spectrum.  However,  the  spectrum  analysis 
of  the  theoretical  seismogram  mentioned  above  gives  a  satisfactory  result,  and  the 
following  things  can  be  learned  from  this  analysis. 


90* 


HOMOGENEOUS  MANTLE 

|T, 


.t# 


iTn 


iTi*  iTif 


iT» 


(1^10) 


Fio.  i. — Spectrum  computed  from  theoretical  seismop 
component  of  displacement  due  to  the  tonional  oscfll 
line*,  the  line  spectra  of  the  free  torsional  oscillation,  ah* 
and  the  intensity  of  the  spectrum  used  for  computinf  d 
Length  of  the  seismograms  used  in  die  analysis  is,  front 
4-0, 7-s.  a-o,  2-0  and  8-0,  respectively,  in  terms  of  the  tint 
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Fig.  i. — Spectrum  computed  from  theoretical  seismograms  (aaimuthal 
component  of  displacement  due  to  die  torsional  oscillation).  Vertical 
lines,  the  line  spectra  of  the  free  torsional  oscillation,  show  the  frequency 
and  the  intensity  of  the  spectrum  used  for  computing  the  seismograms. 
Length  of  the  seismograms  used  in  die  analysis  is,  from  top  to  bottom, 
4-0, 7'2,  a-o,  2-0  snd  8-o,  respectively,  in  terms  of  the  time  unit  (2«s/F,). 
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(1)  Both  frequency  and  the  amplitudes  are  given  accurately.  (See  Figure  1. 
Vertical  lines  are  the  line  spectra  of  the  free  oscillations.) 

(2)  The  resolving  power,  or  the  sharpness  of  the  peaks  depends  on  the  length 
of  the  data  analysed,  which  fact  agrees  with  theory.  However,  the  very  long 
duration  of  data  usually  used*  seems  not  to  be  necessary  for  the  detection  of  large 
peaks,  unless  a  very  good  resolving  sharpness  is  needed. 

(3)  When  there  is  a  large  peak,  the  small  peaks  nearby  are  affected,  and  the 
frequency  and  the  amplitude  are  not  necessarily  given  accurately. 

(4)  The  relative  heights  of  the  spectral  peaks  differ  from  place  to  place.  For 
instance  6  =  60°  and  120°,  which  are  complementary  angles,  have  similar  peaks. 

(5)  The  series  of  peaks  naturally  makes  a  guiding  curve  (see  the  fourth  curve 
in  Figure  1),  which  gives  a  clue  for  the  determination  of  the  mechanism  of  the 
origin.  This  will  be  the  topic  of  a  future  study  (Sat6  &  Usami  1963). 

3.  Group  and  phase  velocities 

By  the  ordinary  method  of  measuring  the  peaks  and  troughs  of  the  observed 
disturbance,  the  apparent  group  velocity  was  obtained  from  the  theoretical  seis¬ 
mogram  (Figure  2).  There  is  no  appreciable  discrepancy  between  the  curve  thus 
obtained  and  the  theoretically  calculated  values  from  the  formula 


U  =  dzjdn.  (3.1) 

In  this  expression  z  is  the  non-dimensional  frequency  of  free  torsional  oscillation 
of  a  homogeneous  elastic  sphere,  which  is  tabulated  in  a  previous  work  of  the 
present  authors  (Sat6  &  Usami  1962). 

Since  the  phase  velocity  and  the  group  velocity  are  connected  by  a  first  order 
differential  equation,  the  phase  velocity  can  be  obtained  by  integrating  the  curve 


Fig.  a. — Graphs  of  group  velocity  (UIV$)  against  period  for  the  case  of 
a  homogeneous  sphere.  Solid  circles  show  the  group  velocity  computed 
by  ordinary  method  of  measuring  die  peaks  and  troughs  of  seismograms. 

Broken  line  is  the  curve  theoretically  obtained  from  (3.1). 

•See  for  example  Benioff  k  others  1961,  Nest  k  others  1961,  Alsop  k  others  1961,  Brunt  k 
others  1961. 
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in  Figure  a.  The  initial  value  at  the  atarting  point,  which  it  unknown,  it  assumed 
to  be  ro8,  109  and  t*io  at  (period)  -  0*04.  All  of  these  curvet  are  shown  in 
Figure  3,  from  which  it  it  possible  to  estimate  the  correct  value  of  the  phase 
velocity. 


Fig.  3. — Phase  velocity  (broken  lines)  obtained  by  integrating  the  group 
velocity  curve  given  in  Fig.  2.  Open  and  solid  circles  are  theoretical 
group  and  phase  velocities  for  the  case  of  s  homogeneous  sphere  com¬ 
puted  from  the  non-dimensional  frequencies  of  the  fundamental  mode 
of  the  free  torsional  oscillation.  Wavelength  is  assumed  to  be  equal  to 

(2»«)/(»l  +  J). 


3 a.  Phase  velocity 

The  phase  velocity  can  also  be  calculated  using  another  method  (Satd,  1955, 
1956,  i960  1960a),  namely,  the  Fourier  analysis  of  two  different  wave  groups, 
say  G 1  and  Gy  This  commonly  used  process  was  employed  in  the  present  study. 
Since  there  are  many  possible  combinations  of  wave  groups,  only  a  number  of 
them  are  illustrated.  Figure  4  contains  seismograms  of  wave  groups  used  in  this 
analysis.  The  calculated  dispersion  curve  (Figure  5)  agrees  well  with  the  theore¬ 
tically  obtained  one,  but  discrepancies  appear  when 

(1)  The  period  (wavelength)  becomes  longer 

(2)  The  difference  of  two  epicentral  distances  becomes  smaller. 

The  relation  between  the  wavelength  and  the  difference  of  epicentral  distances 
at  which  the  discrepancy  becomes  noticeable  (about  0*3  per  cent)  is  briefly  given 
as  follows  (Figure  6  and  the  following  table). 

Difference  of  epicentral  distance  (A)  720°  540°  360°  180°  120°  60° 

Wavelength  (A)  65'  550  42°  34°  23°  120 

Order  of  oscillation  (»)  568  10  15  18 

This  relation  is  summarized  by  the  following  empirical  formula: 

A  —  o-2  A  A  <  180° 

A  -  23+0-06  A  A  >  180° 


(3-a «) 
(3-2b) 
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If  it  ia  assumed  that  the  numerical  computation  for  the  phase  angle  of  the 
Fourier  transform  of  the  disturbance  involves  the  error  of  o»i  rad,  the  first 
formula  can  be  theoretically  deduced.  Any  one  of  a  number  of  effects  can  cause 
the  above  mentioned  error,  but  probably  this  is  close  to  the  error  of  a  single 
observation. 


0-90*  G I  (,WM) 

0-2  v  M 
e-00*  G2(,VM 


Homogeneous  sphere 

0-60*61  (swM) 


0-4  0-7  04 


0>9O*G3  (|WM) 


0-60°  GS  (jW«o) 


II  12 


|.'o  II  1-2  14 

0-90°  G4(|W(o)  „  . 


e-l20°GI  (sWw) 

Hn" 

02  03  04 


IS  14 


8*90*G5  (,W,o) 


1-7 


0*I2O°G3(jW*o) 


10  l-l  1-2  14 


i 

l-S  £0  2-1  2-2 

Fio.  4. — Seismogram*  of  wave  trains  used  in  the  analysis. 


3  b.  Spectrum  and  phase  angle  of  G  waves 

Examples  of  the  spectrum  and  the  phase  angle  of  the  Fourier  transform  of  G 
waves  are  shown  in  Figures  ya  and  7 h.  Both  are  restricted  to  the  phase  Gi  only, 
the  range  of  time  analysed  being  about  0-25,  or  25  min.  The  original  disturbance 
is  shown  in  Figure  4. 

From  the  result  of  analysis  it  is  seen  that: 

(1)  The  spectra  at  two  different  points  are  similar.  (Uppermost  curve  in 
Figures  7 a  and  7 b). 

(2)  There  is  a  zero  point  in  the  spectrum  curve,  which  exactly  coincides  with 
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Figs.  5(a)  to  j(/). — Phase  velocity  (F/K.)  obtained  theoretically  (broken 
line)  from  the  frequency  of  free  torsional  oscillation  and  that  obtained 
by  Fourier  analysis  of  two  different  wave  groups  mentioned  in  each  figure. 

the  zero  point  of  the  Fourier  conjugates  of  the  source  function,  namely 

(sinpfi)/p  =  0,  h  =  0-015.  (3.3) 

(3)  The  distribution  of  the  phase  angle  is  regular,  making  a  smooth  and  almost 
straight  line.  (Middle  curve  in  Figures  7 a  and  76.) 

(4)  The  inclination  of  this  line  is  proportional  to  the  epicentral  distance  as 
predicted  by  the  relation  between  phase  angle,  phase  velocity  and  epicentral 
distance. 

(j)  There  is  a  jump  of  n  on  this  line  at  the  zero  point  of  the  spectrum. 

(6)  Since  we  know  the  phase  velocity,  the  initial  phase  j8 (p)  at  the  origin  is 
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calculated  from  the  following  formula  (Sat6,  loc  cit) ; 

m+m+*rg{f*(p‘,0))-p8IV(p)  (3.4) 

in  which  L  (0)  is  the  polar  phase  lag,  0  the  epicentral  distance  and  V  (p)  is  the  phase 
velocity  as  a  function  of  frequency  p.  The  curve  thus  obtained  also  has  a  jump  of 
it  at  the  zero  point  of  the  spectrum  (Figures  7a  and  7 b  below).  On  the  other  hand 
P(p)  can  be  calculated  as  the  argument  of  the  function  in  (3*3),  which  should  be 
either  zero  (for  the  range  of  the  function  (sin  ph)lP  positive)  or  it  (for  the  range 
of  the  same  function  negative). 


Fio.  6. — Relation  between  the  wavelength  and  the  difference  of  epicen¬ 
tral  distance  at  which  the  error  of  phase  velocity  calculated  by  Fourier 
analysis  becomes  noticeable  (about  0-3  per  cent). 


(7)  Since  there  is  such  a  jump  of  phase  angle,  care  must  be  taken  when  the 
phase  velocity  is  calculated  from  a  single  observation. 

4.  Equivalent  surface  source 

In  the  previous  sections  the  origin  of  disturbance  was  assumed  to  be  on  the 
surface  of  the  sphere.  But  this  is  not  the  case  for  actual  earthquakes  and  the 
question  naturally  arises  what  will  happen  if  the  stress  is  applied  at  a  surface  at 
some  depth  within  the  Earth. 

In  order  to  obtain  the  disturbance  when  the  identical  stress  is  given  on  the 
surface  r  =  b  instead  of  r  —  a,  the  reciprocal  theorem  of  the  field  is  conveniently 
used.  Since  it  was  found  in  the  previous  study  (Sat6,  Usami  &  Ewing  1962)  that 
most  of  the  disturbance  belonging  to  the  surface  wave  is  contributed  by  the 
fundamental  mode,  the  higher  radial  modes  will  be  neglected  and  only  the  funda¬ 
mental  radial  mode  will  be  considered.  From  the  reciprocal  theorem,  it  is  known 
that  the  displacement  at  r  -  b  due  to  the  surface  source 

ip  -  constant  (4.1) 
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Figs.  7(0)  to  7(6). — Spectrum,  phase  angle  of  the  Fourier  transform  of  G 
waves  and  0  (p)  the  initial  phase  at  the  origin. 

is  nothing  but  the  surface  displacement  due  to  a  source  at  a  point  r  =  b.  In  the 
present  case,  however,  the  area  on  which  force  is  applied  is  proportional  to  the 
square  of  the  radial  distance  of  the  source.  Therefore,  (A/a)*  should  be  multiplied 
by  (4.1)  in  order  to  obtain  the  surface  displacement  due  to  a  source  at  r  —  b, 

to  =  constant  '  f*(P )•  (4.2) 

Multiplying  the  factor 

"  la/  /«+i(*a) 


(4-3) 
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by  the  surface  displacement  due  to  a  surface  source  (Sat6,  Usatni  and  Ewing  1962) 

to,  -  constant  •  -f*(p)  (4.4) 

a* 

leads  to  the  expression  (4.2).  This  fact  implies  that  the  surface  displacement  due 
to  a  force  applied  at  r  -  b  having  spectrum  /*  (<)  is  equal  to  the  surface  displace¬ 
ment  due  to  a  surface  force  with  spectrum  /*  (/) .  Bn.  Therefore  Bn  is  termed 
the  “equivalent  surface  source”  which  is  given  in  Figure  8. 


FlC.  8. — Factor  B».  Spectrum  of  the  “equivalent  surface  aource”  is 
obtained  by  multiplying  this  function  by  the  spectrum  of  a  aource  at 
the  depth  (a  — b). 


The  time  distribution  of  the  stress  at  the  origin  was  synthesized  from  the 
spectrum  (sin  ph)/p  and  Bn  .  (sin  ph)lp.  The  first  one  should  give  the  rectangular 
form  of  the  original  source  function.  Because  of  the  neglected  high  frequency 
(maximum  frequency  used  is  63*66),  however,  the  calculated  form  of  the  function 
is  a  little  different  from  the  original  one  (Figures  9 a  and  9 b).  When  the  origin  is 
deep  within  the  Earth,  short  period  surface  waves  have  small  amplitude,  as  is 
easily  expected. 

Summary 

Theoretical  seismograms  obtained  in  previous  studies  of  the  propagation  of 
torsional  disturbance  on  the  surface  of  a  sphere  have  been  used  to  test  and  evaluate 
the  usual  methods  of  analysis  of  seismograms.  Special  attention  is  paid  to  the 
accuracy  of  analysis  and  the  effect  of  parameters  such  as  length  of  record  and 
epicentral  distance.  Both  the  frequencies  and  amplitudes  of  the  component 
spectral  lines  have  been  recovered  by  Fourier  analysis.  Group  and  phase  velocities 
obtained  by  the  peak  and  trough  method  and  by  Fourier  analysis  of  successive 
wave  trains  show  good  agreement  when  compared  with  theoretical  values.  The 
initial  phase,  important  in  detection  problems,  has  been  recovered,  and  agrees 
well  with  theoretical  values.  The  idea  of  an  equivalent  surface  aource  is  discussed, 
using  reciprocity  to  modify  the  spectrum  of  a  source  at  depth. 
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FlO.  9(a). — Time  distribution  of  the  (tret*  at  the  origin  for  the  surface 

source. 


Fic.  9(6).—' Time  distribution  of  the  stress  at  the  origin  for  the  “equivalent 
surface  source,”  which  is  given  by  the  expression  B.  ./•  (p)  (cf.  (4-3)). 
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